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Abstract 

We investigate the hadronic cooling of a quark droplet within a microscopic 
model. The color flux tube approach is used to describe the hadronization of the 
quark phase. The model reproduces experimental particle ratios equally well 
compared to a static thermal hadronic source. Furthermore, the dynamics of 
the decomposition of a quark-gluon plasma is investigated and time dependent 
particle ratios are found. 



1 Introduction 

When the possible existence of a quark-gluon plasma was first postulated in the 
mid seventies ||^ several groups started to investigate unique signatures which 
would allow the identification of such a gas of asymptotically free quarks and 
gluons in heavy ion collisions. Tightly connected to this task was the study of 
the mechanism for the decomposition of a once formed partonic state of matter, 
back to the hadronic world measured in detectors. One was immediately lead 
to the question of the order and signatures of a QGP-hadron phase transition. 
Under the widely used assumption of a first order transition a sharp surface of 
the quark phase could be defined at which the hadronization should take place. 
Applying Gibbs conditions at the phase boundary then gives a macroscopic view 
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on the properties of the initial quark and the final hadron phase. However, the 
dynamical aspects of the plasma disassembly (e. g. the dissociation rate) could 
not be tackled in this simple picture. 

Therefore, several different models for the dynamics of the transition were 
proposed which all focused on the radiation of hadrons from the plasma surface. 
One of these contributions was driven by our late friend Judah Eisenberg 
who proposed a quark-gluon bag which was surrounded by a cloud of pions. 
Both phases interact at the phase boundary via the coupling of the pseudo- 
scalar pion field to the quark vector-field in the interior of the bag. However, 
this and other similar approaches as e.g. § still lacked a dynamical description 
of the conversion process and therefore seemed to be restricted to stationary 
problems in equilibrium. This problem was partially attacked by the applica- 
tion of the color flux tube model to the radiation scenario as proposed in Q . 
Here, the idea of confinement is adopted to describe hadronization dynamics. 
Now a dynamical non-equilibrium calculation of the pion production rates was 
possible. However, all those approaches predicted that the contribution of sur- 
face radiation to the total hadron yield in heavy ion collisions should be small 
compared to the strong hydrodynamic expansion of a hot fireball. Even the 
extension of the flux tube model to four quark flavors and the complete hadron 
zoo could not circumvent this result. 

Nevertheless, the hadron radiation process is yet not completely negligible 
which is underlined by the rough accordance of all those different models within 
one order of magnitude (flg. |l|). But, evidently, the dynamics of the phase 
transition cannot be separated from the dynamics of the QGP phase itself 
which includes the radiation from the surface as well as the expansion. The 
latter, however, requires a new, dynamical deflnition of the "surface" of a QGP 
as this surface will now be time dependent. In the microscopic picture of 
moving quarks (as the flux tube model) a surface cannot be defined consistently: 
it has to be introduced "by hand" as a boundary condition. The complete 
microscopic description will therefore not contain such a surface. Inspired by 
Judah Eisenberg we extended the flux tube approach by assuming the confining 
force acting on an escaping quark as dynamically induced by all the other quarks 
which define the plasma phase, i. e. one quark moves in the effective color field 
of all the others. 

This paper presents this model to describe the dynamics of the hadroniza- 
tion transition. It explicitly takes into account a non-equilibrium evaporation 
scenario. 
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2 Hadron radiation from quark-gluon plasma 



Two fundamental scenarios have been suggested for the decay of a quark-gluon 
plasma within the past fifteen years. The first scenario assumes a hydrody- 
namical expansion of the plasma and a cooling which leads to hadronization 
1^. In a second approach hadronization is constrained to the plasma surface 
mediated by a particular interaction between the two phases. Several hypothe- 
ses have been developed, all of which assume an equilibrated plasma state that 
evaporates mesons according to specific mechanisms. Reciprocally, the plasma 
should also be able to absorb hadrons from outside with a certain probability. 

The pion emission rate of a thermalized quark gas has been calculated [|| by 
introducing a fixed threshold quark momentum for the hadronization of escap- 
ing quarks to pions. Thus the pion yield is determined by the characteristics of 
the initial plasma state, and the energy flux of pions per unit time and surface 
element d^E/d^S dt obeys a Stefan-Boltzmann T'^ law. 

A more sophisticated mechanism of meson evaporation from a plasma state 
uses the cloudy-bag model postulating the existence of a pseudo scalar pion 
field outside the plasma region coupled to the quarks at the surface Q . In this 
model pions are produced by the processes q qn, q ^ qw and qq it ai 
the surface of the bag. The result of this model differs drastically from the 
predictions of |^ , revealing an approximate power law for the emissivity of . 

In this paper we apply a hadronization scheme proposed by Banerjee et 
al. [Q. A quark penetrating the bag surface is pulled back by a strong color 
interaction with the volume of the plasma. The confining force increases with 
the distance from the plasma so that a quark cannot escape. This attractive 
potential is simulated by a color flux tube connecting the emitted quark to the 
plasma by a uniform color field. This tube can decay by creation of virtual qq 
pairs inside the color field. The color charges of the pair produce screening of 
the color field between the initially emitted quark and the quark droplet and 
allow q and q to build a color singlet state (a meson). The remaining q will be 
trapped in a color flux tube again. 

Baryon production can be implemented e. g. by inclusion of (anti-)diquark 
states (uu, ud, . . .) which later may combine with a single quark to form a 
baryon. However, the probability for qq — qq production is strongly suppressed 
by the relatively large masses of the diquark states, so that the predicted baryon 
multiplicities underestimate the experimental result. Therefore, our calculation 
treats diquarks as constituents of the quark-gluon plasma. This mechanism will 
drastically increase the number of emitted baryons. 

The existence of bound diquarks states in the QGP conflicts with the ex- 
pectations of an asymptotically free quark-gluon plasma. However, recent cal- 
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culations of several groups hint at the existence of correlations in the QGP even 
for T<2Tc . 

We assume that a quark-gluon plasma bag has a radius of i? ~ 2 — 6 fm 
and is in complete thermal and chemical equilibrium at temperatures T 
100 — 250 MeV. It is filled with asymptotically free quarks and antiquarks. In 
the following we neglect the gluonic degrees of freedom. 

Let us consider a thermal quark plasma bag with a sharp surface. This 
state will then radiate hadrons from the surface which change the energy and 
flavor content of the plasma. We consider that near the surface the particles 
are homogeneously distributed in space, while their momentum distribution is 

approximated by the Boltzmann distribution exp(— y^p^ + m^/T). We assume 

to have current quarks (and diquarks) inside the plasma. In this approximation 
the mean particle number for each species q inside the bag is 

N,{T,V,^^) = g.e'^^^ ^mlTK,[^) , (1) 

where q runs over all possible (anti-)quark and (anti-)diquark states, and gq is 
the degeneracy factor which counts all possible color and spin states. 

In contrast to the other models of hadron radiation cited in this article, we 
are, in principle, not restricted to a thermal initialization of the plasma phase. 

2.1 Hadronization via fission of color fiux tubes 

The string mechanism proposed in |^ describes the fragmentation of a color 
flux tube created between two point-like quarks excited by a high impact p + p 
(or heavy ion) collision. It is assumed that longitudinal momenta of quarks 
are much higher than the transverse components so that a 1 + 1 dimensional 
description is appropriate. A string stretched between a quark and a large 
bag is different since the color field is not necessarily constrained to a tube. 
However, assuming free motion of color charges on the surface of the bag, we 
assume a flux tube which is perpendicular to the bag surface and will follow a 
possible quark motion perpendicular to the tube. 

A quark crossing the surface of the bag travels outward in accordance with 
the attractive potential of the color flux tube. In 1 + 1 dimensions this potential 
is taken to be linear 

V{z) = K-z, (2) 

which assumes a constant chromoelectric field strength density dV/dz = k, k, 
0.9 GeV/fm. This value can be extracted from fits to Regge trajectories or from 
lattice calculations 0. 
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The total probability of vacuum decay is then given by 



7rm f n 

f=q,qq n=l 

where the sum runs over all quarks and diquarks, p/ is the probability per unit 
four- volume to produce a virtual qfQf pair in a constant field ^, ^ , and my is 
the mass of the produced quarks (or diquarks) . The total probability of vacuum 
decay is then given by This result has the same form as Schwinger's formula 
for e+e~ production in an constant external field Nevertheless, it should 
be noted that the above result is valid only for a static, homogeneous color 
field, i. e. for a flux tube with infinite transverse dimensions. Finite flux tubes 
impose surface effects which may noticeably modify the field configuration and 
pair production probabilities. In |Q it was shown that a string radius of 0.5 fm 
reduces the pair creation rate by about 40%. 

For the quark masses we assume the following values of m„ — rrid = 
5MeV, ms — 150MeV, TOc — 1.5 GeV. For the mass of the diquarks two dif- 
ferent possibilities are taken into account: masses derived from fits to experi- 
mental e+e~ data which allow to extract a probability ratio p{qq)/p{q) = 0.065 
for diquark production and, secondly, diquark masses calculated within 
the Nambu- Jona-Lasinio model [|l3| . One obtains different masses for different 
flavor and spin combinations. Combining m, d and s flavors and using simple 
symmetry arguments for the diquark wave functions (see e.g.Q) one obtains 
three antisymmetric diquark states of spin 

{ud)o — —={ud — du), {us)o — —={us — su), {ds)o — —j={ds — sd) , (4) 
y2 v2 v2 

and six symmetric combinations 

{uu)i — uu, {ud)i ^ —={ud + du), {dd)i = dd, (5) 

v2 

(ss)i = ss, {us)i — —={us + su), {ds)i — —={ds + sd) , (6) 

V 2 V 2 

with spin 1. Since only the anti-triplet representation of the SU{3)c x SU{3)c 
may contribute to a diquark state, we obtain a color degeneracy of 3. Consid- 
ering all the spin and isospin combinations, we arrive at the total degeneracy 
factors for each kind of diquarks given in table 0. 
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(w)o,(ot)o 




{qs)o,{qs)o 


{qs)i,{qs)i 


(ss)i,(ss)i 


LUND mass (MeV) 


420 


490 


590 


640 


790 


NJL mass (MeV) 


234 


822 


537 


962 


1087 


degeneracy 


3 


27 


6 


18 


9 



Table 1: Diquark degeneracies and masses. Masses are taken from fits to the LUND 
string model and from calculations within the NJL model. 

Assuming a string cross section A and a constant break-up probability p in- 
side the string volume we obtain a time-dependent probability for string fission 

t 

-pA J dt' z{t') 

P{t) = 1 - e « (7) 

where z{t) is the distance of the quark from the surface. The flavor of the cre- 
ated quarks and therefore the species of the hadron to be produced is chosen 
according to the probability given by (||), defining all quantum numbers. If this 
choice proves to be energetically impossible the produced flavor will be rese- 
lected. The remaining quark is propagated within the attractive color potential 
from the position of the point of break-up until the flux tube breaks up again 
or the quark is pulled back into the plasma. String fission is only possible while 
the quark is moving outward (p^ > 0) . 

The flavor and the quantum numbers of the produced qq (or qq — qq) pair 
determine the species of the produced hadron. Quantum mechanical bind- 
ing properties have been neglected as well as possible mass reduction of the 
hadrons due to the dense surrounding medium. Thus the free vacuum mass for 
the hadrons is assumed. The probability for the qq pair to have a transverse 
momentum p± is 

nEj 

P{p^) = e ~ with = .y/Af2 + (pJ^)2 . (8) 

From momentum conservation during the string fission we obtain directly the 
transverse momentum of the meson. The longitudinal hadron momentum is 
determined by using the fragmentation function / according to [p^ ^ ^ 

z/(z)^(l-z)"e-^^l/% (9) 

where z = Phadion/Pinitiaiquark ^ud a, /3 are free parameters which are fixed in 
the LUND model to a = 1 and /? = 0.7 GeV"^ 
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Figure 1 : Pion surface bright- 
ness of a baryon-free QGP 
as a function of temperature. 
Three upper lines are the re- 
sults from Danos and Rafelski 
1^, the dotted line are calcu- 
lated by Miiller and Eisenberg 
Q while the three lower lines 
represents the results from 
Banerjee et al. |^] for vari- 
ous model parameters. The 
filled squares are the calcu- 
lations with inclusion of all 
hadrons species while the open 
squares give the results, when 
only pion production is taken 
into account. 



2.2 Results 

Figure |l| shows the surface brightness (^Et^I^S dt of the pion radiation as func- 
tion of temperature T. We compare predictions of three different models by 
Danos and Rafelski (upper lines) , Miiller and Eisenberg |^ (dotted line) and 
Banerjee et al. I (lower lines) with our results. We predict the same orders of 
magnitude as in [ ^ . The bending of the curve increases heavily with increasing 
temperature which is due to the opening of new hadronic degrees of freedom 
at higher energies. To see this, we recalculated the distribution neglecting all 
hadrons other than pions (□) which then gives a good agreement to the former 
results. 

As describes so far, our model assumes a static plasma source. It has been 
shown in ref. ||l7|| that hadron ratios change with time, i. e. during the reaction 
some particle species are produced earlier than others. Furthermore, the hadron 
evaporation in each time step will influence the statistical properties of both 
the plasma phase and the hadronic phase. This dynamical evolution will pull 
the system from an initial fs = into a domain of finite strangeness in both 
phases [ p7| . 

It is interesting that the model is able to reproduce experimental data fairly 



7 



10' 

S 10 



# experiment 

O microscopic (LUND) 

X microscopic (NJL) 

□ Spieles et al. 



i 6 



hadron species 



Figure 2: Final particle ratios in the 
non-equilibrium scenario with initial 
conditions A^g'' = 100, S/A'"" = 
45, /i"" = 0, B^^^ = 235 MeV 
corresponding to the calculations in 
p?! (squares). The circles represent 
the results assuming LUND diquark 
masses jl6j, while the crosses de- 
note calculations for diquark masses 
extracted from the NJL model 
Data from various experiments as 
compiled in |18| are also shown. 



good even in a static approach, as can be seen in fig. g. The initial conditions 
of the plasma phase have been chosen to be equivalent to those used in 
A^" = 100, 5/^'"" = 45, /™* = 0, 5^/4 = 235 MeV. The NJL diquark 
masses seem to reproduce the experimental data better than LUND masses. 
This is not surprising, since the LUND masses have been extracted from fits 
of the string model to high energy data {^/s « 4GeV) while average quark 
momenta in our scenario are certainly below 1 GeV. 

One sees that our calculations fit the data almost as well as the dynamical 
model and are in accordance to a static thermal fit. This is ensured by the 
inclusion of diquarks in our plasma. First, they are responsible for creating 
baryon numbers but secondly they reduce the distill effect of net strangeness. 
Indeed, at finite an excess of q compared to q leads to enhancement of s and 
therefore in the hadronic phase. The qq and qs diquarks also exceed the 
number of qq and (suppressed) diquarks, but the latter may now combine 
with a s to form a A, S or S. This then leads to an enhancement of s in 
the hadronic phase. The net strangeness enrichment in the hadronic sector is 
therefore balanced by the diquark contribution which is mostly influenced by 
the initial baryon chemical potential /Xg. 

This is illustrated in fig. ^ where the strangeness fraction fs — (Ns — N-s)/A 
is plotted as a function of the chemical potential /ig for /is = (diamonds) 
and /is = 25 MeV (circles) . For small diquark production in the plasma is 
suppressed due to the high diquark masses and therefore s is enriched on the 
hadronic phase giving rise to negative fs- For higher /iq the strangeness fraction 
goes to zero which indicates that both effects almost cancel. At /ig — > there 
should be /s = due to restoration of baryon number symmetry. From the 
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Figure 3: (Left) Strangeness fraction fs = [N^ — N-) /A as a function of initial baryon 
chemical potential /i, from microscopic calculation for /i'/"* = (O) and /i'/"* = 
25MeV (O)- (R-ight) Time evolution of chemical potentials fiq and fia assuming 
a totally equilibrated plasma phase at every stage of the reaction with an initial 
temperature T — 190 MeV. The dashed lines give the respective /i's for a plasma 
without diquarks. 



microscopic calculation it is not clear whether there is a smooth transition in 
fs from finite to /iq = or a discontinuity as net strangeness and baryon 
numbers converge to zero. The large deviation of two points at small /Zq indicate 
a significant uncertainty in this observable. 

For small but finite initial /Zs we observe a positive fs for small /Zq, as it 
can be expected. For higher /iq the described s transport mechanism becomes 
more important and finally dominates the initial strangeness excess. 

From fig. I (left) it is obvious that one feature of a dynamical treatment 
of hadronization cannot be described in a static scenario: the accumulation of 
anti-strangeness in the hadronic phase causes an excess of strangeness inside the 
plasma which will then change the strange chemical potential to finite fig > 0. 
The energy stuffed into strangeness will then cause a strong reduction of /x^. 
Figure ^ (left) shows that for fiq > 50 MeV this will cause an even higher 
strangeness enrichment. For small /Zq, on the other hand, the growing /i^ will 
finally pull the system in a domain of positive fs , which means s-dominance in 
the hadronic phase which can no longer be compensated by the small baryon 
chemical potential. 

Fig. ^ (right) elucidates this scenario. The time evolution of fj,s and Hq is 
plotted for initial conditions T = 190 MeV, = 200 MeV and = for 
both a plasma with (solid lines) and without (dotted lines) initial diquarks. 
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Figure 4: Hadron ratios as a function of fiq for = (O) and Hs ~ 25MeV (O)- 
The grey bars indicate the experimental data including error bars. 

The time dependence is obtained by adjusting the thermodynamical conditions 
inside the plasma to the average baryon number, strangeness and energy trans- 
port through the plasma surface calculated by the microscopic model for given 
initial conditions. The assumption of a complete statistical equilibrium of the 
plasma after each time step allows to extract new parameters fj,q, iis (as we are 
treating only particle rates the volume V of the plasma is of no importance for 
particle ratios, unlike the absolute number of produced hadrons which is pro- 
portional to the plasma surface). The temperature is assumed to be constant. 
It is worth to note that in this scenario we do not need to require equilibration 
for the hadronic phase or validity of Gibbs conditions at the transition. Thus, 
our model is a non-equilibrium description of the hadronization procedure. 

In absence of diquarks in the plasma the strange chemical potential mono- 
tonically rises with time. This is in agreement to the results of ref. If 
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diquarks are taken into account the picture changes drastically. After a short 
period the contribution saturates at about 25 MeV. On the contrary the 
fjLq distribution decreases linearly in time which gives possibility to use /iq as 
a time evolution parameter. The strangeness distillery is strongly reduced by 
including diquarks. 

Finally, fig. ^ shows the jiq dependence of different particle ratios compared 
with experimental data for /i^ = and /i^ = 25 MeV, where the latter represents 
the saturation value for /ig in the time dependence. Within the statistical 
fluctuations of the ratios, the two curves cannot be distinguished. Therefore, 
the distributions may be viewed as some kind of time evolution of the particle 
ratios. 

3 The quark molecular dynamics (qMD) 

A novel extension of the quark dynamics can be realized by treating quarks as 
semi-classical particles which interact according to a two-body color potential 
(for details see [Q). This potential is phenomenologically motivated in order 
to mimic the soft gluonic part of a quark-gluon plasma. The Hamiltonian reads 



W = E 2 E ^S-n|r. - r,|) 

where N is the number of quarks. Here we include four quark flavors (u, d, s, c) 
with current masses to„ = = 10 MeV, = 150 MeV and nic = 1.5 GeV. 

Cj?- indicates the color factor which regulates the sign and relative strength 
of the interaction between any two quarks and antiquarks depending on the 
color combination of each pair. They can easily be obtained from the quark- 
gluon vertex factors in color space and read 

Cc \ ^ \a \a 

a=3,8 

in the abelian approximation. This yields e. g. 

^HM = ^1 ' C'flG = +2 ' ^RR = +-'^ ' ^RG = ^2 ■ 

It is worth to note that the relative strength of the above color factors is 
rigorously enforced by the requirement of color neutrality of widely separated 
qq and qqq states. 
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To provide confining properties V{r) is taken to be linear at large distances. 
At short distances the strong coupling constant as becomes small compared to 
1 which causes the one-gluon exchange terms to dominate the interaction and 
therefore induces a coulomb- like potential. In total, one obtains the well known 
Cornell-potential |Q 

, , 3 tts 

V[r) — \- Kr . 

4 r 

In the static case of infinite quark masses this interquark potential has been 
confirmed by lattice calculations over a wide range of quark distances [Q . 

The second request to the model is to define a criterion how to map those 
bound quark states to hadrons. Such a mechanism is essential as the Hamilto- 
nian is not tuned to describe bound and truly confined hadron states. This has 
been done by the requirement that the total color interaction from a pair (or a 
three particle state) of quarks with the remaining system vanishes. Then, these 
qq- and qqq-states do no longer contribute to the color interaction of the quark 
gas. If a bound quark state fulfills the hadronization criterion it will be mapped 
to an appropriate hadronic state with identical quantum numbers. Spin and 
isospin of the hadron is randomly chosen according to the probabilities given 
by the Clebsch-Gordon coefficients while the mass of the produced hadron is 
determined by energy and momentum conservation. 

3.1 Dynamics of the model 

We can now study the dynamical evolution of the system. We use the Metropo- 
lis algorithm to generate a thermal distribution of quarks where longitudinally, 
a Bjorken-like velocity profile Vz = z/t has been imposed. 

Within about 12 fm/c most of the initial quarks are completely hadronized. 
It becomes obvious that this hadronization process is not compatible to a ra- 
diating hadron source, but is ruptured into several small quark blobs before 
finally dissolving. 

As the geometry of the quark phase resembles those of an assumed fire- 
ball in a S-|-Au collision, we compare in fig. || the calculated hadron yields at 
mid-rapidity to those measured in CERN-SPS experiments (see compilation in 
|l8|| ). We find a very good agreement in all MM, MB and BB ratios, while 
the antibaryons seem to be clearly under-predicted. Particle ratios, however, 
proved not to be a very sensitive observable to test the quality of theoretical 
models. Fits of a pure hadron gas |l^ proved to describe data with a compara- 
ble precision as other thermal or hydrodynamical approaches including a QGP 
phase transition or several microscopic simulations as UrQMD ||l^. However, 
the analysis of event-by-event fiuctuations pO| and of the dynamical properties 
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Figure 5: Final state hadron 
ratios from qMD calculations 
(open circles) compared to 
S+Au data at 200AGeV (full 



hadron species circles, taken from ]lq]). 



of the system may yield new insight. Thus, we shall further scrutinize the dy- 
namical properties of the qMD approach. As the initial conditions were taken 
from the idealized Bjorkcn scenario we shall focus to the transverse dimension 
at mid-rapidity. Fig. ^ (right) depicts the transverse mass spectra of various 
hadron species as produced by qMD with the above initial conditions. The dis- 
tribution of the initialized quarks (left part) shows T = 200 MeV for all u, d and 
s-quarks and underlines the initialization procedure. In the final curves, how- 
ever, we observe a significant dependence of the slope of the distribution from 
the mass of the regarded hadron species. We stress again that no transversal 
flow had been imposed. This result is in strong contrast to the fireball model 
proposed in which explicitly required a transverse velocity profile in order 
to obtain such transversal behavior. 



4 Conclusion 

Experimental particle ratios at SPS energies had been shown to be compatible 
with the scenario of of a static, thermally and chemically equilibrated hadron 
gas. The present paper demonstrates that the experimental data are also com- 
patible with a scenario with a quark gluon plasma consecutively evaporating 
hadrons, provided that the plasma contains thermally distributed bound di- 
quark states. Without diquarks the model fails to describe the observed baryon 
abundances. Neither a bag constant nor equilibration of the hadronic phase is 
assumed in this model. 

Under the assumption that the plasma phase itself remains equilibrated dur- 
ing the reaction, a quasi time evolution of the plasma blob can be determined. 
The model predicts strangeness distillation in the quark-hadron transition. The 
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Figure 6: (Left) Transverse mass spectra of the initial thermalized u, d (circles) and s- 
quarks (open diamonds). Both distributions show T — 200 MeV as initialized. (Right) 
Transverse mass spectra of final tt, K, p and A in arbitrary units. One obtains rising 
inverse slope parameters with increasing particles masses. 

inclusion of diquarks reduces the strangeness enhancement. 

The extension of the model to non-equilibrium initial conditions results in 
a proper description of the transverse and longitudinal expansion with a good 
agreement of the hadro-chemical properties. 
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